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ABSTRACT 


In  this  report  we  present  two  alternate  derivations  of  the 
equations  governing  the  thermodynamics  for  a  simple  model  of  fireball 
rise  and  entrainment.  These  two  derivations  place  emphasis  on  different 
sets  of  thermodynamic  variables.  It  is  shown  that,  except  for  difference 
in  notation,  the  resulting  equations  governing  the  fireball  thermo¬ 
dynamics  arc  identical.  One  of  the  reasons  for  presenting  the  two 
different  derivations  is  to  clarify  some  of  the  differences  and  simi¬ 
larities  of  work  published  by  various  authors  in  this  area.  Unlike  the 
work  most  recently  published  by  us,  there  arc  some  other  treatments 
where  terms  involving  the  time  rate  of  change  of  the  specific  heat  ratio 
arc  neglected,  he  point  out  that  these  terms  arc  as  important  as  any  of 
the  other  terms  in  the  equations  governing  the  fireball  thermodynamics 
in  the  cases  where  (1)  energy  storage  in  internal  degrees  of  freedom 
is  significant  and  (2)  where  energy  transfer  via  dissociation  and 
association  processes  is  significant. 
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1. 


INTRODUCTION 


In  a  previous  report* we  presented  results  of  a  computer  code 
designed  to  calculate  the  rise,  entrainment,  and  chemistry  of  a  nuclear 
fireball.  In  that  report  we  also  presented  derivations  of,  (1)  the 
equations  of  notion  of  the  fireball  rise,  (2)  the  equations  governing 
the  thermodynamics  and  (3)  the  equations  governing  the  chemistry  of  the 
fireball.  Kc  also  presented  a  description  of  the  computer  code  U3cd 
to  integrate  the  coupled  set  of  differential  equations  governing  the 
rise,  thermodynamics  and  chemistry  of  the  fireball  during  its  rise 
through  the  atmosphere. 

In  this  report  we  shall  rc-dcrivc  the  equations  governing 
the  thermodynamics  of  the  fireball  using  a  different  notation.  The 
notation  used  here  is  complementary  to  the  notation  used  in  Reference  1 
in  the  sense  that  the  two  notations  tend  to  place  emphasis  on  different 
physical  quantititics  in  somewhat  the  same  way  as  do  the  different  sets 
of  thermodynamic  variables  in  ordinary  thermodynamics ,  which  sets  arc 
related  to  one  another  by  various  Legendre  transformations. 

In  Section  2  we  review  the  derivation  of  the  thermodynamic 
equations  given  in  Reference  1,  in  order  to  make  clear  the  premises  on 
which  that  derivation  was  based.  In  Section  3  we  present  an  analogous 
derivation  based  on  another  thermodynamic  notation  which  is  that  often 
used  in  the  description  of  high  temperature  air  in  equilibrium.  This 

-  J.>»,  5 

type  of  notation  har  been  used  by  various  authors  in  fireball  calcula- 

* 

tions.  In  Section  4  we  shall  show  that  the  equations  derived  in 
Sections  2  and  3  arc  identical  except  for  notation.  In  Section  4  we 
also  present  some  concluding  remarks. 


* 

References  3,  4,  and  5  arc  meant  to  be  typical,  not  exhaustive. 
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2. 


REVIEW  OF  THE  DERIVATION  OF  EQUATIONS  GOVERNING  THE 
THERMODYNAMICS  OF  THE  FIREBALL 


In  Reference  1  we  presented  a  derivation  of  the  equations 
governing  the  thermodynamics  of  the  fireball.  This  derivation  was 
based  on  several  premises  regarding  the  behavior  of  the  fireball  gas 
during  entrninnent  of  cold  air  from  the  surrounding  ambient  atmosphere. 

The  first  premise  is  that  the  fireball  gas  is  <n  kinetic 
equilibrium  although  not,  in  general,  in  chemical  equilibrium. 

The  second  premise  is  that  the  pressure  contributions  of 
the  various  chemical  species  is  additive,  i.e.,  that  Dalton’s  law  of 
partial  pressures  is  valid. 

The  third  premise  is  that  ».hc  contributions  to  the  interna) 
energy  of  the  fireball  gas  from  each  of  the  chemical  species  is 
additive. 


The  fourth  premise  is  that  the  equations  governing  the 
thermodynamics  of  the  fireball  during  rise  and  entrainment,  can  be 
derived  by  consideration  of  two  infinitesimal  processes.  The  first 
infinitesimal  process  is  the  mixing  of  an  increment  of  ambient  air 
with  the  fireball  at  constant  pressure.  The  second  infinitesimal 
process  is  an  adiabatic  expansion  of  the  fireball. 

The  first,  second,  and  third  premises  allow  us  to  express 
the  pressure  and  internal  energy  (refered  to  unit  volume)  of  the 
fireball  gas  and  the  ambient  atmosphere  in  the  form 

(2.1) 
(2.2) 


P  *  Z  n.kT  ,  E f  «  Z  n.[f.(T)  ♦  C  ] 
i  l  r  i  1  1  1 

P  -  InjkTa  «  Ea  *  +  Cjl 
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In  tqs . (2. 1 )  and  (2.2)  the  sum  over  i  is  n  sum  crcr  the  wny  chemical 
species  of  the  firchall  while  the  sum  over  j  is  a  sum  over  the  few 
chemical  sjiccies  of  the  surrounding  ambient  atmosphere.  denotes 

the  number  density  (concentration)  of  the  i**1  species,  T  the  fire¬ 
ball  (kinetic)  temperature,  T  the  temperature  of  the  ambient 

A 

atmosphere,  k  holtrnann's  constant,  f  ^  (T)  the  temperature  dependence 
of  the  internal  energy  contribution  of  the  i1^  species  and  the 

heat  of  formation  of  the  i**1  species. 

In  Lqs.(2.1)  and  (2.2)  the  pressures  of  the  fireball  and  the 
ambient  atmosphere  arc  the  same,  corresponding  to  pressure  equilibrium 
of  the  fireball  with  the  ambient  atmosphere. 

In  Inference  1  we  combined  the  two  infinitesimal  processes 
mentioned  above  and  derived  the  equations  which  govern  the  thermo¬ 
dynamics  of  the  fireball.  To  avoid  any  confusion  that  might  arise 
in  the  reader's  mind  about  the  effects  arising  from  each  of  the 
processes  we  shall  here  consider  each  process  separately  and  show  that 
the  same  results  arc  obtained. 


Kc  first  consider  the  process  of  mixing  at  constant  pressure 
a  mass  AM  of  ambient  air  with  the  fireball  gas.  Kc  can  associate  a 
volume  AVfl  of  the  ambient  air  parcel  and  we  have 


(2.3) 


where 


Vy.n . 

i  J  J 


(2.4) 
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In  Iq.('.-l)  represents 

molecular  weight  of  the  ;l ^ 
process  we  have 


the  density  of  ambient  air  and  the 

species.  For  the  constant  pressure 


1  final 
1  initial 


Vfinal 


(V  ♦  A  V  | )  £  (n4  •  An^JlfjO  •  ATj)  .  CJ  (2.  Sj 

V  »'l  I f i  (T)  *  S1  *  A\*>jlfJ(V  *  V 

C.fc) 

V  ♦  AV,  (2.7\ 


'initial 


V  ♦  AV 

a 


C.«) 


In  l.qs,  (J.S)  -  (2,5),  f  denotes  the  total  internal  energy  of  the 
fireball,  AVt  the  change  in  fireball  volume,  An^  the  change  in 
species  concentrations,  and  AT,  the  change  in  the  fireball  tempera¬ 
ture.  The  constant  pressure  process  is  one  in  which  no  heat  is  added 
tc  the  thermodynamic  system  consisting  of  the  fireball  and  the  ambient 
air  parcel  of  mass  AM.  fhus,  the  first  law  of  thermodynamics  yields 


0 


dQ  * 


Vfinal 


initial  ^  *  ‘final  *  initial 


(2.9) 


he  substitute  the  expressions  (2.S)  -  (2.8)  into  l.q.(2.9)  and  keep 
through  first  order  infinitesimals.  The  result  is 


AVilp  *  I!  n4  (f4*Ct)  J  ♦  V  £  An^f^Cj)  ♦  V  £  -jf  AT, 

AM 

•  ^  IP  •  “  nj(fj(T,)  *  Cj)l 

where  f.  and  dfj/dT  are  evaluated  at  T.  The  equation  of  state 
yields 

Ap,  -  0  -  £  AnJkT  ♦  £  n  kAT, 
i  1  i  1 


(2.10) 


(2.11) 
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1 


Eqs.  (2.10)  and  (2.11)  contain  the  information  obtained  by  the  consi¬ 
deration  of  the  infinitesimal  process  of  mixing  at  constant  pressure. 


We  now  consider  the  second  process  of  adiabatic  expansion 
of  the  fireball.  We  have 

Efinal  =  (V  *  AV*  *  Av*)  E  (n.  ♦  An.1  ♦  An.*)lf.(T  ♦  AT,  ♦  AT2)  ♦  C.] 

(2.12) 

(2.13) 

(2.14) 

(2.15) 


initial  "  (V  4  AV»)  Z  ("i  ♦  An‘.)(f.(T  ♦AT,)  ♦  C.) 
Vfinal  "  V  ♦  AV,  ♦  AV2 


initial  *  V  *  AV‘ 


where  the  subscript  2  and  superscript  2  refer  to  the  adiabatic  expansion 
process.  The  first  law  of  thermodynamics  becomes,  for  this  adiabatic  process 

0  -  dQ  -  P(Vfinal  -  vinitial)  ♦  Efina|  *  ^initial  (2.1<>) 

Substituting  l:qs.(2.12)  -  (2.15)  into  Eq.(2.1b)  and  maki  ng  use  of 
Lq.(2.10),  we  obtain 

df . 


AVj[p  ♦  £  n.(f.  ♦  C.)l  ♦  V  £  An.2  (f.  ♦  C.)  -  V  £  n.  AT2 


(2.17) 


The  equation  of  state  yields 


P  ♦  Ap2  *  Z  (nA  ♦  An^  ♦  An^)  k  (T  ♦  AT,  ♦  AT2) 
or  with  Eq. (2.11) 

Ap2  ■  Z  An.’kT  ♦  £n.kAT2 
i  1  i  1 


(2.18) 


(2.19) 


Eqs. (2. 17)  and  (2.19)  constitute  the  information  provided  by  the 
adiabatic  expansion  process. 


L 
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Adding  i'.qs.(2.II)  ,md  (2.10)  and  noting 


Ani  •  Ani*  ♦  An,2  ,  AT  ■  AT \  *  AT  2 

wc  obtain  the  differential  for*  of 

P  ■  £  n.kT 
i  1 


(2.20) 


(2.21) 


i.c.,  the  equation  of  state.  Adding  liqs.(2.10)  and  (2.17)  and  noting 


AV  ■  AVj  ♦  AVj  , 


the  liq.  (2. 4)  and 


M  ■  V  £  P.n. 

i  1  1 


(2.22) 


(2.23) 


wc  obtain 


o  ■  Jglp*  piffi-S))  -Cj) 

)) 


dn. 


df 


i  dT 


'Slu.-A 

1  dM /  i  1  1  ' 


Hr, L1 — J  tP  ♦  rn  <f  .c 


(2.24) 


Eq.(2.24)  is  identical  with  Eq.(3-45)  of  Reference  1 


Lq . (2.21)  and  (2.24),  which  govern  the  thcnsodynanics  of  the 
fireball,  arc  identical  with  those  obtained  by  coobining  the  two  in¬ 
finitesimal  processes  (as  was  done  in  the  derivation  given  in  Reference 
1).  Since  coobining  the  processes  shortens  the  derivation  wc  shall 
treat  the  problem  that  way  in  the  following  section. 
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3. 


ALTERNATE  DERIVATION  OF  THE  THERMODYNAMICAL  EQUATIONS 


In  this  section  we  shall  derive  the  equations  governing  the 
thermodynamics  of  the  fireball,  using  a  different  notation  than  that 
of  the  previous  section.  In  this  notation  we  write  the  equation  of 
state  in  the  font; 


pV  * 


(3.1) 


where  Z  is  the  number  of  moles  of  air  at  temperature  T  per  mole 
of  air  at  STP  at  sea  level.  V  is  the  volume  of  the  fireball,  M 
the  mass  of  the  fireball,  T  the  temperature  of  the  fireball,  R 
the  universal  gas  constant,  mo  the  gram-inolccular-mass  of  STP  air 
at  sea  level,  and  N  the  number  of  moles  in  the  fireball.  l:q.(3.1) 
is  an  equation  of  state  for  a  real  (as  distinct  from  an  ideal)  gas  and 
therefore  the  quantities  Z  and  N  arc  not  constant  but  vary  when  dis¬ 
sociation  and  association  processes  occur.  The  gas  constant  is  given  by 

R  =  AK  (3 ,2] 

where  A  is  Avogadro's  number  and  k  is  Boltzmann's  constant.  The 
gram-molccular-mass  is  given  by 

m0  =  A  p\  (3.3) 

where  ik  is  the  average  molecular  mass  associated  with  STP  air  at 
sea  level. 


Ke  write  the  internal  energy  of  the  fireball  in  the  form 


H  =  NCy(T)T 

where  the  number  of  moles  of  gas  in  the  fireball  is  given  by 

N  ■  2  (£)  ■  vz  w.> 

The  enthalpy,  ||,  of  the  fireball  is  written 
II  «  H  ♦  pV  =  NCp(T)T 


(3.4) 


(3.5) 


(3.6) 
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Note  that  the  quantities  Cv  and  Cp  of  llqs. (3.4)  and  (3.6)  are 
not  the  specific  heats  per  mole  as  they  would  he  conventionally  defined 
in  thermodynamics.  In  the  co'.vcnt ionnl  sense, 


C 


V 


(3.7) 


where  Cv',  c;  arc  the  specific  heats  per  mole?  Kc  shall  find  it 
convenient  in  what  follows  to  define  a  quantity 


Y  =  Cp(T)/Cv(T) 


(3. 8j 


where  C  Cv  arc  defined  by  Eqs.(3.4)  and 
specific  heat  ratio  y' ,  on  the  other  hand, 


(3.6).  T’ c  conventional 
is  given  by  y'  =  C,'/Cv'. 


he  now  consider  the  combined  infinitesimal  process  of  adding, 
at  constant  pressure,  a  parcel  of  ambient  air  of  mass,  AM,  to  the  fireball 
and  then  expanding  the  fireball  adiabatically.  From  hq . (3.1)  we  obtain 

■  “(£)TP  *  Z(£H  *  z(i)TAp  (3.9) 

The  initial  internal  energy  of  the  system  consisting  of  the  fireball  and 
the  ambient  air  parcel  is 


E.  ...  , 
initial 


NCV(T)T  .  4NaCv(Ti)T: 


(3.10) 


where 


AN 

a 


„  /AM\ 

Z  =  Z  AV  (p  /m0) 
a\m0  f  a  a  a 


(3.11) 
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The  final  internal  energy  is 


‘■final  =  (N  *  AN)  CV  (T  *  ATHT  +  AT) 


The  initial  and  final  volumes  of  the  system  arc 


V.  .  .  , 
initial 

^filial 


V  ♦  AV 

a 

V  ♦  AV 


(3.12) 


(3.13) 

(3.14) 


The  first  law  of  thermodynamics  for  this  adiabatic  process  is 


0  =  dQ 


P^Vfinal  "  'initial^  +  ,;final 


E. 


initial 


(3. IS) 


he  substitute  Kqs.  (3. 10) ,  (3.12),  (3.13)  and  (3.14)  into  Eq.(3.1S)  and 
keep  through  first  order  in  infinitesimals.  The  result  is 


rdC* 

p(AV  -  A\'n)  ♦  ANCyT  ♦  N,  dT 


T  -  C. 


AT  *  “aWb 


(3.  ICO 


Kq . (3.9)  can  be  rewritten  (using  Hq . (3. 1)) 

ApV  ♦  pAV  =  ANRT  ♦  NRAT  (3.17) 

Multiplying  hq.(3.17)  by  C^(T)/R  and  subtracting  the  result  from 
Eq. (3. 16)  we  obtain 

CV  dC 

p(AV  -  A\'a)  ♦  —  (ApV  ♦  pAV)  ♦  N  — -  TAT  =  ANaCV(Ta)Ta  (3.18) 


Ke  now  introduce  the  quantity 

C  ♦  R 

Y  =  - - 

V 

and  note  that  since  p  =  M/V,  we  also  have 

AV  _  AM  Ap 
V  =  M  ‘  p 


(3.19) 


(3.20) 
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Kith  these  relations,  liq.  (3.11)  and  the  relation 


Y 


a 


W  ♦  R 

C(Ta) 


(3.21) 


* 

we  find  that  liq.  (3. 18)  becomes 

1  dO  „  1  I  dj>  ^  1  dM  /  y  -  1  VM  ZaTa)  1  dy  1 

p  dt  y  p  dt  M  dt  |  \  Y;i  -  1  /\  Y  /  ZT  j  y  dt  (y  -  1) 


(3.22) 


This  equation  and  hq.  (3.1),  which  we  write  in  the  form, 


P  = 


(3.23) 


arc  the  analogues  of  liqs.(2.24)  and  (2.21)  of  the  previous  section. 
liqs.(3.22)  and  (3.25)  govern  the  thermodynamic  behavior  of  the  fireball. 


Alternate  forms  of  Eqs.(3.22)  and  (3.23)  can  be  used  instead 
but  such  alternate  sets  of  two  equations  for  the  thermodynamic  behavior 
must  be  equivalent  to  liqs.(3.22)  and  (3.23). 


Note  the  presence  of  the  last  term  on  the  right  hand  side  of 
Bq.(3.22)  involving  the  time  derivative  of  the  quantity  y. 

As  we  shall  see  in  the  following  section,  this  term  is  of  the 
same  size  as  the  other  terms  appearing  in  Kq.(3.22). 
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4. 


EQUIVALENCE  OF  THE  EQUATIONS  GOVERNING  THE  THERMODYNAMICAL 
BEHAVIOR  OF  THE  FIREBALL 


In  this  section  we  shall  show  that  the  set  of  i.qs.  (-.21) 
and  (2.24)  is  identical,  except  for  notation,  with  the  set  {.iven  by 
Fqs.(3.22)  and  (3.23).  For  this  purpose  we  find  it  convenient  to  re 
introduce  the  variables  M,  V  in  place  of  p  in  Hq.(3.22).  From 
Eqs.(3.20)  and  (5.22)  we  obtain 


0 


1  dV 
V  dt  + 


1_  l_  d£  l_  dM 
Y  P  dt  M  dt 


Z  T 
a  a 

ZT 


I  I 

Y  dt  (y-1) 


(4.1) 


l:rom  the  Fqs.(2.21),  (3.1)  and  (3.2)  we  obtain 


V  £  n. 


(4.2) 


From  Eqs .(2.1)  and  (3.4) 


SCVT  =  vEn^f.  (T)  .  c.) 


(4.3) 


and  substituting  fo_*  N  from  Kq.(4.2)  we  get 


?  "i  (fi  *  c0 

CV(T)  ■  *  ly- 


(4.4) 


From  Eq .(3.1)  and  the  analogous  equation  for  the  ambient  atmosphere 


p  =  Z  i— )  pT  =  Z  (— )  p  T 
1  \mo'  r  a  \mo/  a  s 


(4.5) 
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or 


Z  T 
a 

ZT 


a 


£uini 

YTT. 

j  J  J 


(4.6) 


The  last  form  of  Eq.(4.o)  follows  from  the  relation  (2.4)  and  an 
analogous  relation  for  p. 


It  is  convenient  to  use  the  notation  of  Eqs.(2.1)  and  (2.2), 
and  in  the  following  we  drop  the  subscript  f  from  the  internal  energy 
per  unit  volume,  E^.,  of  the  fireball.  The  equations  of  state  (2.21) 
and  (3.1)  arc  obviously  identical  if  we  mal:c  the  identification  (4.2). 
Using  Hqs.(4.2),  (3.3)  and  (3.5)  we  obtain 


Z 


77  y  n. 
V3  i  ] 

*>ini 


(4.7) 


It  remains  to  show  that  Eqs.(2.24)  and  (4.1)  arc  identical.  Kith  the 
notation  of  Eqs.(2.1)  and  (2.2)  we  can  write  Eq.(2.25)  in  the  form 


0 


s 


I  dv 

V  dt 


(P  ♦  r.) 


I  dM  _P_, 
M  dt  P„  U 

U 


(4.8) 


From  Eq.(4.4)  we  have 


Thus 


cvO)  =  Ak  ^  ,  Cp(T)  =  Cy  ♦  R  =  Ak  (1  ♦E/p)  (4.9) 


Y 


E_L 1 

i:  > 


p  ♦  E 
r  a 


E 

a 


(4.10) 
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and  therefore 


vi  „  Y  -1 

Y- 1  _  p  a 


p-h:  *  y. 


_1L_ 

p*i; 


(4.11) 


hith  this  notation  l.cj .(4.1)  becomes 


0  «  V  .  (A.)  i  &  .  I  dM  j>  l 

V  dt  \P+1:/  P  M  dt  P„  \  P+L 


p.l  \  *  I:  -  i!i  p 

—i)  .  l!_  <lt  (4.12) 


p  (p*ro 


Die  second  term  on  the  right-hand-side  is  cancelled  by  the  first  member 
of  the  last  term  on  the  right-hand-side  and  we  obtain 


1  df. 


v  eTt  +  (p*r)  dt 


1  <?M  p  /P*Ea\ 
m  dt  pa  y  P+i  j 


(4.13) 


which  is  identical  with  Hq.(4.S).  This  concludes  the  demonstration  of 
equivalence. 


finally  we  remark  that  the  last  term  on  the  right-hand-side 
of  hq .(4.12)  is,  in  general,  of  the  same  magnitude  as  the  other  terms 
in  Kq .  (4.12).  flic  last  term  arises  from  the  term  involving  dy/dt  in 
i-q.H.l).  This  term  is  evidently  neglected  in  some  other  treatments  of 
the  fireball  thermodynamics.  The  last  term  in  Kq . (4.12)  vanishes  in 
the  case  of  an  ideal  gas,  where 


_  _  NRI  r  1  T 

P  v  *  F:  *  v  NLVT 


(4.14) 


and  the  quantities  N  and  Cy  arc  constant.  In  this  ease  we  have 


d£  B  NR  £T  NRT  1_  dV 

dt  "  V  dt  V  V  dt 


(4.15) 
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cP 

i!t 


(4.16) 


Ni'V  ill  NCV  1  dV 

V  ilt  "  V  i  dt 


so  th.it 


ilil  .  i!il 

dt  dt 


P 


0 


(4 


The  treatment  presented  in  the  previous  sections  assumes  that  the  fire¬ 
ball  and  ambient  atmosphere  gases  have  two  properties  in  common  with  an 
ideal  gas.  Ihcsc  properties  are  (1)  that  the  contributions  to  the 
Pressure  from  the  various  species  are  additive  and  (2)  that  the  contri¬ 
butions  to  the  internal  energy  from  the  various  species  are  additive. 
1-ror  th  point  of  view  of  a  microscopic  picture  embodied  in  a  statisti¬ 
cal  mechanics  description,  these  assumptions  about  the  bulk  properties 
of  the  gas  arc  tantamount  to  the  assumptions  (1)  that  the  behavior  of 
the  gas  can  be  characterized  in  terms  of  a  temperature  and  (2)  that 
at  constant  temperature  and  species  concentration  (i.e.,  when  there  are 
no  chemical  reactions)  the  interactions  among  unbound  particles 
(potential  energy  associated  with  interactions)  make  a  negligible  con¬ 
tribution  to  the  internal  energy  compared  with  the  energy  associated 
with  the  kinetic  and  internal  degrees  of  freedom  of  the  particles  of 
which  the  gas  is  composed. 


On  the  other  hand,  the  treatment  assumes  (A)  that  the  con¬ 
tribution  to  the  internal  energy  from  the  internal  degrees  of  freedom 
of  the  various  species  is  important  and  (B)  that  the  change  in  the 
internal  energy  of  the  gas  due  to  changes  in  the  species  concentrations 
is  an  important  effect.  The  ki.ictic  energy  released  when  recombination 
takes  place  and  the  kinetic  energy  absorbed  when  dissociation  occurs 
provide  important  energy  transfer  mechanisms.  Suffice  it  to  say  that 
the  last  term  on  the  right-hand  side  of  Eq.  (4.12)  makes  an  important 
contribution  to  both  effects  (A)  and  (B). 


•  17) 
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